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Abstract 

A new kind of phase transition is proposed for lattice fermion systems with simplified f 2 config- 
urations at each site. The free energy of the model is computed in the mean-field approximation for 
both the itinerant state with the Kondo screening, and a localized state with the crystalline electric 
field (CEF) singlet at each site. The presence of a first-order phase transition is demonstrated in 
which the itinerant state changes into the localized state toward lower temperatures. In the half-filled 
case, the insulating state at high temperatures changes into a metallic state, in marked contrast with 
the Mott transition in the Hubbard model. For comparison, corresponding states are discussed for 
the two-impurity Kondo system with f 1 configuration at each site. 

1 Introduction 

In some uranium compounds with 5/ 2 configuration (U 4+ ) the CEF ground-state can be a nonmagnetic 
singlet. The CEF singlet is also realized in some praseodymium compounds with 4/ 2 configuration 
(Pr 3+ ). In these cases the spin entropy of the system can go to zero as temperature decreases even 
though interactions with conduction electrons or with f electrons at other sites are absent. This is in 
striking contrast with the case of cerium compounds with 4/ 1 configuration (Ce 3+ ); the entropy does 
not disappear at zero temperature if a Ce ion is isolated because of the Kramers degeneracy associated 
with the f 1 configuration. As a result the system chooses, depending on the interaction between Ce sites, 
among a magnetically ordered state, a Fermi liquid state, a superconducting state, and so on in which 
the entropy vanishes at zero temperature. 

In these lattice fermion systems, which we call the f 2 lattice hereafter, the itinerant state is also 
possible if the hybridization is large enough. Thus both the localized f-electron picture and the band 
picture can be a starting point to understand the actual compounds with f 2 configuration. The most 
interesting situation occurs when the energy scale of the CEF singlet state is comparable to that of the 
itinerant state. Then both states compete for the stability. 

Suppose we have the CEF singlet as the ground state of the f 2 lattice, but its energy is only a 
little lower than the itinerant state. If the itinerant state is metallic, the entropy increases linearly as 
temperature increases. The temperature scale here is the Kondo temperature, and is related to the large 
density of states at the Fermi surface. On the other hand increase of the f-electron part of the entropy 
in the CEF state follows the exponential law, and is much less significant in a low temperature range. 
Thus there is a possibility for a phase transition to occur from the CEF singlet state to the itinerant 
state as temperature increases in the f 2 lattice system. Even if the itinerant state is a Kondo insulator, 
the entropy can increase more rapidly than that in the CEF state since the energy gap decreases with 
temperature. In the latter case the system changes from a metal to an insulator as temperature increases. 
This is opposite to the case of the Mott transition where the low-temperature phase is an insulator. 

The purpose of this paper is to demonstrate the presence of a phase transition between the CEF singlet 
and itinerant states in the f 2 lattice system at zero and finite temperatures. As the first step to explore 
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an f 2 lattice system, we take the simplest possible approach and apply the mean field approximation with 
account of both itinerant and localized characters of f electrons. The plan of the paper is as follows: In 
the next section, we introduce the model and derive the mean-field equations. The same approximation 
scheme is applied to the two-impurity Kondo system with f 1 configuration at each site in Section 3. It 
turns out helpful to compare the electronic states of both models. The relative stability of different phases 
in both the f 2 lattice system and the two-impurity Kondo system is studied in Section 4. Presence of 
the phase transition at finite temperatures is demonstrated. The final section is devoted to discussion of 
results with attention to possible experimental relevance. 



2 Model and Mean-Field Equations 

We introduce an f 2 lattice model as follows: 

ker v=oc,{3 
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where i is the site index, and v and /i are channels of conduction and f electrons, respectively. We 
express the CEF singlet and triplet using a pseudo-spin operator of f electrons for each channel: = 

(1/2) J2-y8 fifij^ysfind: where a is the vector composed of Pauli matrices. The spin operator of conduction 

electrons is given by S?, = (1/2) J2-yS c w^^fS<HvS- ln e q- (0); the second term with J > gives 
antifcrromagnetic interaction between f and conduction electrons on each site. In the presence of the 
potential scattering term (l/4)n^n^, the Kondo scale Tk = D exp(— 1/ Jp c o) is reproduced correctly in 
the mean-field approximation. Here, D is a half width of a conduction band and p C Q is the density of 
states per spin of conduction electrons at the Fermi level. We note that the sum of the spin exchange 
and potential scattering terms is half of the permutation operator whose eigenvalue is 1 for the quasi-spin 
triplet and is —1 for the singlet. The last term with / > in eq. ([!]) represents the CEF splitting. This 
splitting / is also correctly given by the mean-field approximation due to the term (1/4)7% n^. The 
restriction = 1 is imposed on eq. (|l|) to simulate the strong Coulomb repulsion between f electrons. 
We take a mean field as 

for fi = 1,2 and v — a, (3. This mean field represents the fictitious hybridization between f and conduction 
electrons. We say "fictitious hybridization" in the sense that the real hybridization is absent in the model 
with fixed occupation of f states. However wc neglect in this paper the phase fluctuation which makes the 
mean field vanish. The physical motivation for the neglect will be discussed in the final section. Another 
mean-field is given by 

Ri = ^ y E(fLfaa) (3) 

(T 

which expresses the mixing between two f orbitals on each site. This mixing gives rise to bonding- 
antibonding splitting of localized levels. In eqs. @ and (||) we assume that the mixing is allowed only 
for the same spin directions. The Lagrange multiplier terms 



E E 

=1,2 



are added to eq. (]l|) to enforce the constraints on the number of f electrons. In the mean-field approxi- 
mation the number of f electron per site and channel is fixed only as average. Therefore care is necessary 
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about spurious charge fluctuations included. We discuss this aspect of the mean-field theory again in the 
final section of the paper. 

Assuming equivalence of different sites and channels, we put e^ v = £k, Vi^ u = V^„ = \V\ exp(— iO^), 
Ri = R and Ai M = A. Setting the origin of energy at the Fermi level, we write the free energy per site as 



r EE trln 
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where N s is the total number of the sites, if = J + 7/4 — A and A = irp c0 \V | 2 . The 2x2 matrix M with 
Alatsubara frequency iu n has components 

2\V\ 2 

M MM = -i Wn + i ( + ^-LJ — {n = 1,2), (5) 
iui n - e k 

2\V\ 2 

M 12 = -R+ , IM 5, (6) 

VjJn - £k 
2\V\ 2 

M 2 i = -R* + - ' 1 B*, (7) 

VjJn - £k 

where S = {exp(— i<f> a ) + exp(— i<j>p)}/2 with <p v = 6\ v — 9 2l j (y = a, (3). Here we represent the magnitude 
of the hybridization between different channels of f electrons via conduction electrons by the parameter 
S. Let us take the bases in which the on-site hybridization between f and conduction electrons occurs 
only with the same channel. Even for this case, the intersite hybridization between the channels 1 and 
2 of f electrons can occur. This is because the point group symmetry around each site is not relevant to 
intersite interactions. From a detailed analysis we find that R and S can be chosen real. Then we can 
assume — 1 < S < 1 in the following. 

We derive the mean-field equations by requiring the free energy to be stationary against variation of 
A, if and R. As a result three characteristic states appear: First the Kondo state is the itinerant state 
where f-electrons hybridize with conduction electrons (A ^ 0, R = 0). Secondly the CEF state is the 
localized state where f-clcctrons form the singlet of quasi-spins at each site (A = 0, R ^ 0). Thirdly 
the mixed state has a character interpolating between the Kondo and CEF states (A ^ 0, R ^ 0). 
For the CEF state we always have the solution A = if = and R = I{f(—R) — f(R)}/2. Here, 
f(w) = l/{exp(w/T) + 1} is the Fermi distribution function with energy w. 

In the following we consider the case where the number 7V C of conduction electrons is twice the 
number of lattice sites, and the conduction bands without hybridization have constant density of states 
Pco = l/(2D) between the band edges ±D. Then the system has the insulating ground state if f electrons 
form energy bands, since f and conduction electrons have the hybridization gap at the Fermi level. This 
is called the Kondo insulator. If f electrons are localized, on the other hand, the Fermi level lies in the 
middle of the conduction band, and the system becomes metallic. For the Kondo and mixed states, 
each conduction band is split into two pieces with a band gap between them. We use the notation 
E c = e f - (-l) p(c) J R, and S c = 1 + (-l) p ^Z where p(a) = 1, p(b) = with C = a,b. The band edges 
after splitting are given by 

Die = 

D 2<: = 
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Then the mean-field equations are given by 
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We solve eqs.(]8|)-(|TC]) for various values of dimensionless parameters I/Tk and S at zero and finite 
temperatures. We note that the mixed state is metallic because the condition RA ^ requires a finite 
density of states between bonding and antibonding f levels. 

3 Two-Impurity Kondo System 

In the course of understanding the electronic state in the f 2 lattice system, a necessary step is to clarify 
the difference from the f 2 impurity system. This impurity system is further related to the two-impurity 
Kondo system with f 1 configuration at each site. Namely, the f 2 impurity is considered as the short- 
distance limit of two Kondo impurities. Fortunately we have detailed knowledge about the two-impurity 
system by the mean-field theory [[l], the Quantum Monte Carlo and the numerical renormalization 
group H H . In this paper we derive the ground state and the free energy in the same level of approximation 
as is used for the f 2 lattice system. Then by comparing the electronic state of the f 2 lattice system and 
that of the impurity system, we obtain information about the influence of the lattice periodicity. 

In ref. [Q the two-impurity Anderson model was solved by the mean-field theory. It was shown 
that the intersite hybridization gives smooth change from the limit of two independent Kondo states 
(Kondo pair) to the pair singlet state as the intersite interaction increases. Physically we expect the 
same situation even though the occupation of f electrons at each site is very close to unity. The limiting 
case is described by the two-impurity Kondo model. Although the Anderson lattice model is more general 
than the Kondo lattice model, the CEF state is harder to treat in the mean-field theory. Since we have 
adopted the Kondo lattice model with f 2 configurations, we need to solve the two-impurity Kondo model 
for comparison. 

The two-impurity model is given by 



where J, I > and j(= 1, 2) labels sites of f electrons. There is only a single conduction band since even 
in this case different screening channels are present around each impurity. We take the mean-fields in the 
form analogous to the f 2 lattice system: One is given by 



which represents the fictitious hybridization between f and conduction electrons at each site. The other 
is given by 
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which expresses the mixing between two f electrons. As before the Lagrange multiplier terms 



3=1.2 



are added to eq. ( |TT| ) to enforce the constraints on the number of f electrons. 

There are two dimensionless parameters A, B which represent the intersite hybridization effect via 
conduction electrons. Namely we define 



Here, the magnitudes of A and B depend on both the distance r between f sites and the band structure of 
conduction electrons, but they are always less than unity Q]. The parameter A causes asymmetry in the 
density of states of bonding and antibonding f states: In the case R > 0, with A > (< 0) the density of 
states of the bonding states becomes wider (narrower) than that of the antibonding states. On the other 
hand the parameter B controls the splitting between bonding and antibonding f states. The left-hand 
side of eq. ( p^ ) is analogous to the term with 3 in eq. (^) if one interchanges sites in the former with 
channels in the latter. In the two-impurity system there are three characteristic states: The first is the 
Kondo-pair state where the Kondo effect occurs independently at each site (A ^ 0, R — 0); the second is 
the pair-singlet state where the pair-singlet of f electrons is formed without help of conduction electrons 
(A = 0, R 7^ 0); the third is the mixed state which interpolates the above two states (A ^ 0, R ^ 0). 

4 Stability of Itinerant and Localized States 
4.1 Zero temperature 

We have solved the mean-field equations numerically at zero temperature both for the f 2 lattice and the 
two-impurity systems. Table |l| summarizes the parameters used in the calculation. Figure |l](a) shows 
the ground-state energy per site in the f 2 lattice system at zero temperature. The origin of energy is 
taken to be that of the Fermi sea without f electrons. The abscissa represents the bare CEF splitting 
in units of Tk- The notations Ek, -Ecef and E m i x represent the ground-state energies of Kondo, CEF 
and mixed states, respectively. The effect of intersite hybridization 2 depends only on its absolute value. 
Thus results of E m - lx with 5 = and 3 = 0.4 are shown as representative cases in Fig. |l|(a). 

We find that E m { x is larger than Ek and -Ecef for all combinations of parameters I /Tk and 2. 
Therefore the change from the Kondo state to the CEF state occurs discontinuously at the critical point 
J/Tk = 4. The mixed state which would have interpolated the Kondo and CEF states smoothly is not 
stabilized actually; with increasing intersite hybridization, the mixed state with energy E m { x becomes 
larger in the mean-field theory. This is seen by the fact that i? m i x with 5 = 0.4 is larger than that with 

5 = in Fig. |l|(a). The reason is the following: If 3 > and R > 0, the density of states of the bonding 
f states has a larger width than that of the antibonding f states. The ground-state energy is given by 
the sum of single-particle energies of occupied states. Namely, we integrate the total density of states 
multiplied by w from — oo to 0. Since the integral without w is fixed by the number conservation, the 
total energy increases by the asymmetry induced by 5. Similarly in the case of 3 < and R > 0, the 
asymmetry of the density of states in the opposite direction increases the energy again. 



For comparison, Figjl](b) shows the ground-state energy of f electrons in the two-impurity system at 

zero temperature. We have tried various values of A and confirmed that A does not influence the relative 
stability of the phases. On the contrary the value of B drastically affects the ground state. Hence, we fix 
A = —0.2 and vary B as a free parameter. We note that if B = 0, E m i x is larger than both Ek and £ff pan- 
tor any value of I /Tk- This situation is analogous to that in the f 2 lattice system. As a result an abrupt 
change from the Kondo-pair state to the pair-singlet state occurs as I /Tk is increased. At the critical 
point of I /Tk = 2.5, the two kinds of singlet states are degenerate. Thus one observes the divergence 
of physical quantities such as the susceptibility and the specific heat coefficient. We have checked that 
this level-crossing behavior remains the same as long as \B\ < 1/ir. On the contrary, if the hybridization 
effect is large (l/ir < \B\ < 1), the Kondo-pair state connects continuously with the pair-singlet state 




(12) 
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through the mixed state. In this case no divergence occurs. This is shown in Fig.|l|(b) by the result that 
E m i x (B = 0.4) is lower than both Ek and -Eff pa ir for any value of I/Tk- These results obtained in the 
mean-field approximation agree with those in refs. ||, ||, [|] for the two-impurity Anderson model. 

We note that A controls asymmetry of the density of states of f electrons with respect to the Fermi 
level, just as S does in the f 2 lattice system. We have confirmed that -E m ix becomes larger as \A\ increases 
with B being fixed. Thus the parameter A plays the same role as 5 in the f 2 lattice system. However in 
the lattice system we do not have the parameter corresponding to B in the impurity system. Hence there 
is no stable mixed state in the f 2 lattice system in our calculation. Mathematically the ineffectiveness S 
in stabilizing the mixed state comes from the absence of k-summation in the self-energy of the f 2 lattice 
system. 

4.2 Finite Temperature 

The mean-field equations are solved numerically also at finite temperatures, and the free energies are 
derived. Figure ||(a) shows temperature dependence of free energies per site for three different states in 
the f 2 lattice system: Kondo, CEF and mixed states. Even though the CEF singlet is the ground state, 
there is a case where the itinerant state is realized at higher temperatures. We find that the free energy 
of the mixed state is larger than those of the other two states for all values of parameters I/Tk and 5. 
Therefore, the transition between the Kondo and CEF phases occurs as a first-order one. 

Figure ||(b) shows free energies in the two-impurity system. The notations Fk, ^ffpair and F m i x 
represent the free energies of the Kondo-pair, pair-singlet and mixed states, respectively. It is seen that 
with B = 0, -Fmix is larger than Fk and impair- As in the case of zero temperature, the mixed state is not 
stabilized as long as \B\ < 0.08 for all values of I/Tk- On the other hand, F m i X with B — 0.4 in Fig. ^(b) 
is lower than both Fk and -Fcef- We- have checked that with 0.08 < \B\ < 1 the mixed state is stabilized 
at all temperatures. 

From these results we infer that the parameter B in the two-impurity system plays a decisive role 
also at finite temperatures. Furthermore in the f 2 lattice system there is no temperature region where 
the mixed state is stabilized. In other words, the hybridization effect which mediates between Kondo and 
CEF phases is ineffective, and the transition occurs discontinuously. 



5 Discussions 

5.1 Comparison with two-impurity systems 

In considering the relevance of the mean-field theory, we first take the case of impurity systems. The 
physical difference between the two-impurity Kondo and Anderson models is whether there is charge 
fluctuations of f electrons or not. Reliable knowledge is available for both models from several numerical 
calculations. Computation using the numerical renormalization group derived a level crossing between 
the Kondo-pair state and the pair-singlet state [|| . As a result divergence of the staggered susceptibility 
occurs at zero temperature. On the contrary, a quantum Monte Carlo calculation for the two-impurity 
Anderson model p| found continuous behavior in physical quantities. This apparent conflict was resolved 
by Sakai et al. Q who identified the origin of the continuous crossover as the bonding-antibonding 
splitting of f orbitals. In the Kondo model, the splitting is absent because there is no charge degrees of 
freedom for f electrons. Thus the divergent behavior is purely a formal consequence of the model since 
there should always be some amount of charge fluctuations in real systems. 

For our purpose of studying a new type of phase transition, we regard our model given by eq.(|l|) 
only as a simplified form of Anderson-type models which are more difficult to analyze by the mean-field 
theory. Then the finite order parameters for various phases are rather to be regarded as properties of a 
corresponding Anderson- type model. In formally exact treatment of eq. (|l]) , all of our order parameters 
would vanish identically in contrast with the results of the mean-field theory. However by the same 
reason as explained above for the two-impurity models, we would rather accept the results of the mean- 
field theory as a physically possible consequence for more realistic models. 
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5.2 Comparison with the f 1 lattice system 

It is instructive to take the limits I — ► and S — > in eq.(Q). Then the system becomes equivalent to 
two independent Kondo lattices. In the half-filled case the ground state is either the Kondo insulator or 
a magnetically ordered phase. The latter state is due to the RKKY interaction which is not taken into 
account by the mean-field theory. We note that the ordered state can be either metallic or insulating 
depending on the magnetic structure. If it is ferromagnetic, the half-filled conduction band leads to the 
metallic state. As one increases J from zero, the ground state should change from the magnetically 
ordered state to the insulating one. Subsequently the nature of f electrons changes from the localized 
character to the itinerant one. Since the entropy is different in the two phases, a phase transition can 
occur from one phase to the other as a function of temperature. 

Let us compare this phase transition with another one which is known as an artifact of the mean-field 
theory. Namely, as temperature increases in the mean-field theory, the order parameter A of the Kondo 
insulating phase decreases continuously to zero around the temperature T ~ 7k. In the exact theory the 
local gauge fluctuation washes away the transition completely. We emphasize that the possible transition 
between the Kondo insulator and the magnetically ordered phase should survive the fluctuation effect. 

Now we consider the case of finite / and S. In the mean-field theory, the Kondo insulating phase does 
not feel the effect of I and S. The resultant state is the same as the direct product of the two f 1 lattices. 
Hence the second-order transition around T ~ Tk is again fictitious. One can ask at zero temperature 
how the magnetic state changes as I increases continuously from 0. For small /, we have the ordered 
induced moment which arises by mixing of the singlet and triplet levels. At certain critical value of /, 
the CEF singlet will become more stable than the induced moment. The situation is analogous to the 
spin chain problem where exchange interactions of alternating strength form dimers with an excitation 
gap. We note that the f electrons are always localized for any value of I. Thus the change to the Kondo 
insulator can occur as a phase transition, although both states are spin singlets. 

5.3 Effects of charge fluctuation on the phase transition 

In actual f 2 lattice systems there should always be charge fluctuation as discussed above. Then any 
exact eigenstate has some amount of hybridization between f and conduction electrons. However, there 
can still be two different kinds of hybridized states: The first one can be reached by perturbation theory 
with respect to hybridization. This state is connected with the localized f 2 state. The other state is the 
itinerant state which is not accessible by such perturbation theory. The latter state instead is simply 
described by the band picture of f electrons. Thus possibility of the phase transition remains even though 
effects of charge fluctuations are included. 

We note that the transition to the CEF singlet phase is of first order. A first-order transition should 
be less sensitive to fluctuation effects than a second-order one. We plan to check the robustness of the 
phase transition by using theories ||, ^, |?j more reliable than the mean-field theory. 

5.4 Possible experimental relevance 

With respect to experimental relevance, we have to consider also the case where the number of conduction 
electrons deviates from 2N S . In the Anderson lattice model with dominant f 2 configurations, the itinerant 
state then has a finite density of states of f electrons at the Fermi level, and hence is metallic. In this 
case the Fermi level is shifted from the center of the hybridization gap, and the average occupation of 
f electrons per site also deviates from 2. Thus in reality the transition from the itinerant state to the 
localized one is not always an insulator-metal transition. 

Concerning possible relevance of our theory, we mention two uranium compounds: UNiSn and 
URu2Si2- In the former case, the insulating state at high temperatures changes via a first-order tran- 
sition into a metallic state at T — 43 K. In contrast to our model, however, the metallic state shows 
the antiferromagnetic order ||. As long as the localized picture applies to the low-temperature phase, 
the driving force of the transition may be similar to the one discussed in this paper. It is necessary to 
include the induced moment for more detailed analysis of UNiSn. In the latter case of URu2Si2, the CEF 
singlet model accounts for gross features of highly anisotropic susceptibility and metamagnetic transition 
. The high-temperature phase is metallic showing the Kondo effect in the resistivity. A clear anomaly 
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in the specific heat is observed at temperature To = 17.5 K |L0|, and the resistivity shows the metallic 
behavior also below Tq. By neutron scattering |rT| the antiferromagnetically ordered magnetic moments 
were observed below To. The magnitude of the moment is only 0.04^b which is smaller by two orders 
of magnitude than the usual magnitude observed in similar compounds UT2Si2(T = Pd,Rh). Moreover, 
the growth of moments with decreasing temperature does not follow the mean-field behavior. Strangely, 
the NMR does not probe the internal field below To Thus there is a possibility that the apparent 
antiferromagnetism is not a true long-range order, but due to very slow fluctuation of U moments. 

In any case the specific heat jump at To is too large to be accounted for by the tiny magnetic moment. 
Thus, the proper order parameter in this ordered phase remains to be identified jl3], [u], [if], We note 
that the inelastic neutron scattering probed a feature which looks like a CEF excitation below T . 
This fact may be a key to identify the order parameter. The phase transition seems to be of second order. 

In summary, we have shown that the phase transition from the itinerant state to the CEF singlet state 
occurs as temperature decreases in the / 2 lattice system in the frame of the mean-field theory. Properties 
of the f 2 lattice system were discussed in comparison with the two-impurity system at zero and finite 
temperatures. We suggest that the competition between localized and itinerant states of f electrons is the 
fundamental driving force for phase transitions in some uranium compounds such as UNiSn and URu2Si2- 
It remains to see to what extent the fluctuation effect beyond the mean-field theory affects the phase 
transition. 
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Figure 1: The ground-state energies of (a) the f 2 lattice system, and (b) the two-impurity system. In (a), 
Ek, Eqef and E m { x correspond to the Kondo (itinerant) , CEF (localized) and mixed states, respectively, 
and S is a parameter to characterize the strength of intersite hybridization. In (b), Ek, Sffpair and E m - lx 
show the energies of Kondo-pair, pair-singlet and mixed states, respectively. The parameter B (\B\ < 1) 
characterizes the strength of intersite hybridization effect. Another hybridization parameter A is fixed to 
be —0.2 (see text). 



Figure 2: Free energies of (a) the f 2 lattice system and (b) the two-impurity system. In (a), Fk,Fcef 
and F m i x indicate the Kondo (itinerant), CEF and mixed states, respectively. In (b), Fk, -Fffpair and F m j x 
indicate the Kondo-pair, pair-singlet and mixed states, respectively. The parameter A is fixed to be —0.2 
as in Fig.l. 
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Table 1: Parameters for calculation. 

D T K I 

10 4 [K] 1 [K] ~ 10 [K] -1-1 
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